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Abstract 

A general fusion method to find solutions to the reflection equation in higher 
spin representations starting from the fundamental one is shown. The method 
is illustrated by applying it to obtaining the K diagonal boundary matrices in 
an alternating spin 1/2 and spin 1 chain. The hamiltonian is also given. The 
applicability of the method to higher rank algebras is shown by obtaining the 
K diagonal matrices for a spin chain in the {3*} representation of su(3) from 
the {3} representation. 
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Quantum integrable systems are usually treated by imposing periodic boundary condi- 
tions. Recently there has been a growing interest in exploring other possibilities compatible 
with integrability. There is a method, proposed by Sklyanin in the framework of the al- 



gebraic Bethe ansatz |I[ and relying on previous results by Cherednik [£3] , for obtain 
integrable models with non-periodic boundary conditions. Sklyanin's original formalism, 
which assumes the model is invariant under the parity and time reversal symmetries, was 
extended to more general systems in Refs. [§-[5| among others. A careful analysis within 
the framework of algebraic structures was carried out in Refs. For the XX Z model, 

an alternative approach, based on a diagonalization schema via vertex operators |§, can be 
found in Ref 0. 

Sklyanin's formalism introduces a boundary K + (8) matrix, which must verify the well- 
known reflection condition expressed by 



R{6 - 9') \K + {6) <g> i] R{6 + 6' 



K 



+ (a'\ 



K + (9') <g> /] R(6 + &) \k + {6) <g> i] R(9 - &). 



(1) 



Matrix K + (6) together with R (8) determine the integrable system with open boundary 
conditions. In this way, solutions to models associated to the fundamental representation of 
different algebras have been found in Refs ||Tl~l-|T6| 



The concept of obtaining results in higher dimensional representations from those in lower 
dimensional representations is known as fusion methods. For R matrices it was worked out by 
Karowski |l7j and Kulish, Reshetikhin and Sklyanin |18 and for K matrices by Mezincescu 



and Nepomechie [19| . 

In this paper, we discuss the application of a fusion method to a non-homogeneous 
chain with spin 1/2 and spin 1 in alternating sites. The hamiltonian is also given. As an 
example of extension of the method to higher rank algebras, the method is applied to a chain 
based on the basic representations of sit (3). The notation used is explained in figure 1. 
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Fundamental representation space fs) 
Non-fundamental representation space (a) 





Fig. 1 



In order to find the reflection matrices in a higher representation a, one should solve Eq. 
(H) with all operators in the required representation. This would require the knowledge of 
the ( cr )i?(6 l ) operators, whose complexity increases with the dimension of representation d, 
as they are represented by d 2 x d 2 matrices. Instead, we propose an alternative method. 

Taking into account that Yang-Baxter equations have validity in any representation, 
we can use them with either equal or different representations for the auxiliary space and 
the site space. Thus, we distinguish two representations s and a and consider the transfer 
matrices 

T^f {9) = L<%1 (6) ® LM (0) ® • • • ® (6) , (2) 

and 
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T { :f (9) = L%> (6) ® L<r% (6) ® • • • ® (0) 



(3) 



Following the method in ref. [0, we build the doubled transition operators 



U^\e)=T^(9) {s) KtA9)Tt' 



(4) 



and 



U ( :f (B) = (6) {a) K+ s (6) T^ 1 



(5) 



These expressions are represented in fig. 2. 



(s, s) 

U a , b 0) 




(a. s) 

U a ,p (6) 




Fig. 2 



The [/ operators fulfill the Yang-Baxter equation 



= [f/ (s ' s) (0') ® J]L (CT ' S) (0 + 0') [f/ (,T ' s) (0) ® J]L (s ' ,7) (0 - 0') , 



(6) 



which is graphically expressed by fig. 3, 
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e-e 





Fig. 3 



whereas the reflection matrices K must verify a relation similar to (pi), namely 



(0 _ [ {a)K + (6) ® I]L^ (6 + 0') [ (S) K+ {6') <g> I] 

= [ {S) K + {9') ® I\L™ (6 + 6') [ (a) K + (6) ® (0 - 6') 



(7) 



which is graphically expressed by fig. 4 





Fig. 4 
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The last equation provides the key for finding the matrices ^K + (6) from the matrices 
in the fundamental representation ( s )K + (9) and the operators L^'^ and L^ s ' a \ 

In order to illustrate the method, we apply it to su(2). In this case we start from the 
XX Z model obtained in the fundamental representation of U q (su(2)) and find the reflection 
matrices for spin a. 

The operator 

L^) (0) = sinh (e + -7 (l + 2a 3 ® S 3 )) + sinh 7 (V ® S~ + a~ ® S + ) (8) 



where q = exp7, is already known |TT|. The operator L^a) can be obtained from L^'°^ by 
transposition. The diagonal solution for ^K + {9) can be found in and is given by 

'fsr + sinh(e+ -6) 

K + sinh (e+ + 6 



(9) 



and e + being arbitrary parameters. By introducing (H) and (H) in (0) we obtain 2a 
independent equations that determine the diagonal solution for the ^K + (8) matrix. 



a) K+{6) = S itj K + fl sinh ( e+ - (a + ± - /) 7 + sign(* - ± - Z)fl) 



(10) 



As a particular case, we can take the spin a — 1, then the i^-elements are |20| 

(0) = X + sinh ( e+ - \ - 0) sinh (e+ + | - 0) 
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(i)#£j(0) = ^+ sinh ( e + + ^ - ) sinh ( e + - o 

{1) KU9) = K + sinh ( e+ - J + 0) sinh (e + + 2 + 



(11a) 
(lib) 
(11c) 



2 ' ^ 2 

The hamiltonian can be obtained from the K and ( a )R(6) matrices. For the a = 1 case, it is 

N-l 

H = sinh( 7 ) [SiSi+i ~ {SiS i+1 f + 2sinh 2 ( 7 ) (SfSf +1 - {S*S? +1 ) 2 ) 
i=i 

+2 (1 — cosh( 7 )) ((SfSf + i + Sf Sf +1 )SfS* +1 + S^S^ +1 (S^S^ +1 + 5'fS'| / +1 )^ 

+2sinh 2 ( 7 ) ((^f) 2 + (^) 2 ) ; 

sinh 2 ( 7 ) cosh( 7 ) / . w™\2\ 
+ s 1 nh(e_-|; smh ( e . + i) (-h(2e_)^ - smh( 7 )(^) 2 ) 

SU,h( ~ ' (sinh(2e + )Siv + sinh( 7 )(^) 2 ) + ctn • J (12) 



2 sinh (e+ — |) sinh (e+ + £ 



which corresponds to the hamiltonian given in and |22] with boundary terms [^(J . 

For the isotropic spin 1/2 case (XXX model), the diagonal solution of the reflection 
equations is 

'e+-0 
e + + 9 

Using our procedure, we found the solution for spin a 

2a 

(f , 4- rr + 

2 ° v 2 



. K+(6) 

ISO v y 



w lso ^i(^) = n ( e + + °- + J - fc + ^ - \ - k )°)- 



(13) 



(14) 



i=i 



Of course, this result can also be obtained from [10] by replacing 7 by —77 and by 776* in the 
limit 77 — > 0. 

A very interesting application of the method is to a non-homogeneous chain combining 



different kinds of spin in alternating sites. In reference [2|[] a non-homogeneous chain based 
on the su(2) algebra combining the spin s = 1/2 and a = 1 in alternating sites is described 
. We shall now describe this chain with open boundary conditions. Following the same 



notation as in 23 , we have 



L(|'3)(0) 











f ai 














o\ 


fa 








o\ 
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6_ 
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aij 



(15) 



with 



a{9) =sinh(0 + 7 ), 
b(6) = sinh (6), 
c{6) = sinh (7), 



(16a) 
(16b) 
(16c) 
(16d) 



and 



= sinh(0 + ~ 7 ), (17a) 
b + (6) = sinh (6 + ^), (17b) 
6_(0) = sinh (0-2), (17c) 



ci(0) = sinh (7)^/2 cosh 7. (17d) 
The monodromy matrix, using the representation a = 1/2 as auxiliary space, is 

iT a a J(0, a) = Lli(6) ■ dt 2 {e + <*)••• LtJ^a^WL^^ + a). (18) 
The doubled monodromy matrix can be calculated by using reflection matrix (||) for spin 

1/2 

|E/*(0) = iT a a J(M) • |< d (0) • (19) 

where, for purpose of simplicity, we have made a = 0. Taking into account that 

h K-(p) = tK +t (0-i), (20) 

we can obtain the hamiltonian through the t operator. The explicit form for the hamiltonian 
can be written by dividing it into different parts 

2N-3 2N-1 

H = X! fyi,i+M+2) X! fyt,i+l) + tyl) +^(2AT-l,2AT)- ( 2 1) 

i=l <=1 

iodd iodd 



The first two parts are the usual hamiltonian 2vJ and the others are introduced by the 



boundaries. After a long calculation, for the three-sites operator we find 

h [tJ+i,t+2) = \ ( cosh (7) - + ~ cosh (7) (2 cosh (7) - cosh (3 7 ) - l)a 
+ sinh 2 ( 7 )[cosh cosh(^)U + 4 sinh 2 (2j)U 
+2 sinh (7) sinh (27) W + 2 sinh (-) sinh (27) [cosh (7)] W 
-i sinh 2 (7)0- • S 2 + 4 sinh 2 (2j)a ■ S 2 , (22) 

where 
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a = a i - I ■ a {i+2) + a { ■ I ■ a {i+2) 

U = I-Si +1 a x l+2 ^+I-Sl +1 ( x^ 2 \ 
U = I-Sl +1 a^ +2 \ 

V = 4- {S x , S z }^ ■ a ^ + a y ■ {S y , S z }^ ■ 

+al ■ {S x , S z }^ ■ a x ^ + 4 ■ {S y , S z }^ ■ erf 2 \ 
W = a\- Si m)2 • 4 +2) + at ■ S^ 2 ■ a^ 2) . (23) 

The two-site operator is 



C») = Sinh2 W i COsh ^)[ COsh • + < ■ S v +1) ) 

+ cosh 2 ( 7 )<7* • Si i+1) + sinh 2 ( 7 )/ w • Sf +1)2 ] (24) 



and the two boundary terms 



sinh 2 



^ ) -4^hl 7 T sinh(2 ^ coth(e - ) ^ 1) (25) 



and 



h [m-i,2N) = sinh2 (7) coth (e+) 



2 cosh (7) + cosh (37) {2N+1) 

(J 



4 sinh (7) 



- cosh (7) sinh {i){af N ~ l) ■ Sf N ? + I ■ S^) 

-sinh(|)[cosh( 7 )]M ( xr- 1 ) ■ {S x ,S,}M + ff W-i) . {S y ,S z }^)]. (26) 

In this hamiltonian, the terms proportional to the identity have been removed. 

In the alternating chain, we can build another system by taking the representation a = 1 
as auxiliary space. The doubled monodromy matrix S U is obtained with the reflection matrix 
\K given in ( |lla| -c) and the operators L^ 1 ' 1 ) and L^'a) which are obtained by transposing 
(§) . The new it operator commutes with the previous 1/2^ one, so it follows that both 
hamiltonians also commute. 

As an example of the application of the method to a higher rank algebra, we will find 
the reflection matrices in a spin chain whose site space is in the {3*} representation of su(3), 
starting from the chain whose site space is in the other fundamental representation {3}. 
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We have two solutions in the fundamental representation {3} given in references \12\ and 
I6H. They are, 



{3} 



K+ l (d)= K + e e sinh(e + - 



{3} K+ 2 (9)=K + smh(e + + -6), 
{3} Kl 3 (e)=K + e 2e smh(e + + l- 



(27a) 
(27b) 
(27c) 



and 



{3} 



K+ 1 (6)= K + sinh (e+ 



{3} K+ 2 (6)=K + e 2e smh(e + --6) 
{3} K+ 3 (6)=K + e e smh(e + + h). 



(28a) 
(28b) 
(28c) 



By applying the method described above and knowing that 



L({3},{3*})(0) 
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a J 



(29) 



with 



a{6) 



'3 7 
sinh [ -9 H — 
\2 2 



b{9) = sinh (-(^ + 7) 

9+7 

c(0) = — sinh (7)e _ 2 _ , 
<i(6') = — sinh (7)e~^ _ . 



(30a) 
(30b) 
(30c) 
(30d) 
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we obtain, for the {3*} representation, the diagonal solutions 

m K+ tl {9)= K + e 9 smh(e + + h - |), (31a) 

m K+ 2 (0)= K+e 2 *sinh( e+ -h- 1\ (31b) 

m K+M= K + sinh ( C+ -\e-\) (31c) 

and 

m Kt tl {9)= K+ e 2<> * inh ( e + + \® + \\ (32a) 

m K+ 2 (9)= K + sinh ( C+ + le + 1), (32b) 

{ s-}«JsW= K + e e sinh (e + - ^ + 1). (32c) 

The method is quite general and will be applied in future work to look for non-diagonal 
solutions with arbitrary spin and to models based on higher rank algebras. 
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